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Abstract 

The quantum damped harmonic oscillator is described by the master equation with 
usual Lindblad form. The equation has been solved completely by us in arXiv : 0710.2724 
[quant-ph]. To construct the general solution a few facts of representation theory based 
on the Lie algebra su(l, 1) were used. 

In this paper we treat a general model described by a master equation with generalized 
Lindblad form. Then we examine the algebraic structure related to some Lie algebras and 
construct the interesting approximate solution. 

Quantum Computation (Computer) is one of main subjects in Quantum Physics. To realize 
it we must overcome severe problems arising from Decoherence, so we need to study Quantum 
Open System to control decoherence (if possible). 

This paper is a series of [1] and [2], and we study dynamics of a quantum open system. First 
we explain our purpose in a short manner. See [3] as a general introduction to this subject. 
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We consider a quantum open system S coupled to the environment E. Then the total 
system S + E is described by the Hamiltonian 

H s+E = H s ®l E + ls®H E + H I 

where H$, He are respectively the Hamiltonians of the system and environment, and Hj is the 
Hamiltonian of the interaction. 

Then under several assumptions (see [3]) the reduced dynamics of the system (which is not 
unitary !) is given by the Master Equation 

^p=-i[H s ,p]-V(j>) (1) 

with the dissipator being the usual Lindblad form 

= \ E { A ) A iP + P A ) A i ~ 2 ^M]) • (2) 
{j} 

Here p = p(t) is the density operator (or matrix) of the system. 

Similarly, the equation of quantum damped harmonic oscillator (see [3], Section 3.4.6) is 
given by 

—p = — i[ua^a, p] — — (a^ap + pa) a — 2apa^ — — (aa) p + paa} — 2a) pa)j , (3) 

where a and a^ are the annihilation and creation operators of the system (for example, an 
electro-magnetic field mode in a cavity), and p, v are some real constants depending on the 
system (for example, a damping rate of the cavity mode). 

Since this is one of fundamental equations in quantum theory it is very important to con- 
struct the general solution. In [3] or [1] some methods to construct a solution are presented. 
However, in [2] we gave the general solution in the operator algebra level. This is a very 
important step. 

Our method is as follows : we clarified a certain algebraic structure arising from the Lie 
algebra su(l, 1) and its representation in the equation and constructed the general solution 
by use of the (well-known) disentangling formula. See for example [3] and [B]. The method is 
popular in Quantum Optics, while it may be not in the field of Quantum Open System. 
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In this paper we want to generalize the model in order to examine deeper algebraic struc- 
tures. Namely, we consider a master equation with generalized Lindblad form defined by 



—p = — i[ua)a, p} — ^ (atap + pa} a — 2apa i ) — — (aa) p + paa) — 2a) pa) 
~ (a 2 p + pa 2 - 2apa) - | {{a ] fp + p(a f ) 2 - 2a t pa t ) 



(4) 



where k is a complex constant satisfying the condition pv > \k\ 2 which ensures the positivity. 
See for example [tJ^I . 

Then we examine an algebraic structure related to the Lie algebras sm(1,1) and su(2), and 
construct interesting approximate solutions by use of it. 

In order to solve the equation we use the method in [1] once more. For that we review a 
matrix representation of a and a' on the usual Fock space 



JF = Vect c {|0),|l>,|2>,|3),-.-}; \n) 
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where e ld is some phase. Note that aa) — a) a + 1 = N + 1. 



the paper it is called the Kossakowski-Lindblad form not the generalized Lindblad one. It may be 
suitable. 
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For a matrix X = {x^) G M(T) 



X 



( 

X\\ X\2 ^13 ' 

X21 X22 %23 ' 

^31 ^32 ^33 ' 

v ; ; ; 

we correspond to the vector X e jrdim c T ag 



X — \ x ij) > X — (Xn, X±2, X13, ■ ■ ■ ; X21, X22, X23, ■ ■ ■ 5 ^31) x 32, x 33, 



(7) 



where T means the transpose. The following formula 



AXB = (A® B T )X 



holds for A,B,X e M{T). 

Then flU) is transformed into 



= Hp(t) 



p(t) = e tH p(0) 



(9) 



where 



H = -iuo{N ® 1 - 1 ® iV) 



-^{iV ®l + l®iV-2a® (a f ) T } - ^{(iV + 1) ® 1 + 1 ® (iV + 1) - 2a f ® a T } 
-^{a 2 <8 1 + 1 <g> (a 2 ) T -2a® a T } - ^{(a f ) 2 <8 1 + 1 <8 ((a f ) 2 ) T - 2a f <8 (a t ) T }. (10) 



v 



Moreover it is rewritten as 



~ H-v. . , .N®1 + 1® N + 1®1 t T /■ txT 

if = — — 1 ®l-(p + v) V va) ®a 1 + pa® (a T ) J 



-2iuj- 



N (8 1-1(8 JV _ / t / hr (a f ) 2 ® 1 + 1® ((a f ) 2 ) r 



+k < a (8 a 



a 2 ® 1 + 1 ® (a 2 ) T 



(11) 



Now let us examine the algebraic structure of H. 
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By setting 

K 3 = -(iV®l + l® N+l® 1), K + = a t ®a T , K_ = a ® (a t ) T (12) 

where N T = N, then we have 

[K 3 , K + ] = K+, [K 3 , K_] = -K_, [K + , KJ\ = -2K 3 . (13) 

Namely, {K 3 , K + , K_} is a set of generators of su(l, 1) algebra, [2]. 
By setting 

Js = ^(N® 1-1® N), J+ = a t ®(a t ) T , J_ = a ® a T (14) 

, then we have 

[J 3 ,J + ] = J + , [J 3 ,J_] = -J_ J [J + ,J_] = 2J 3 . (15) 

Namely, { J3, J + , J_} is a set of generators of su(2) algebra. 
By setting 

K 3 = i(JV g> 1 - 1 <8> N), K + = \ {(at) 2 <g> 1 + 1 ® ((a^f } , = i {a 2 ® 1 + 1 ® (a 2 f} 

(16) 

, then we have 

[K 3 , K + ] = K + , [K 3 , K-\ = -K_, [K+, KJ\ = -2K 3 . (17) 

Namely, {K 3 , K + , K_} is a set of generators of su(l, 1) algebra. 
By setting 

L 3 = -(N <8> 1 - 1 <g> AT), L + = J + - K + , L^ = J_- K_ (18) 

, then we have 

[L 3 ,L + ] = L+, [L 3 ,L_]=-L_, [X+,L_]=0. (19) 

We also note that 

[J+,K + ] = [J-,K.] = 0. (20) 
However, {X3, A-} and {L 3 , L + , L_} don't commute except for 

[L 3 ,/?3] = [L3,A' + ] = [L3,A'_]=0, (21) 



see 0. 

A comment is in order. For later convenience let us write down the remaining commutators. 

[Ks, L + ] = -i{(at) 2 ® 1 - 1 ® ((at) 2 ) r }, [K 3 , L_] = I{( a 2 ® 1 - 1 ® (a 2 f}, 
[K+, L+) = a f ® (a f ) T - (a f ) 2 ® 1, [/<+, L_] = a (g) a T - 1 <8> (a 2 ) T , 
L+] = -a f ® (a f ) T + 1 <g> ((a f ) 2 ) T , L_] = -a <g> a T + a 2 ® 1. 

Then f TTTT) is rewritten like 

# = - - U 1 ® 1 - (/x + ^)A' 3 + + \iK_ - 2iuL 3 + kL + + kL_. (22) 

What we want to do is to calculate the evolution operator e tH , which is in general not easy. 
Since {K%, K + , KJ\ and {L3, L+, L_} don't commute it is reasonable to assume 

e tS PS e i£ ^* e *{-(^+i')K- i +^++^} e i{-2^L3+SL + +KL-} (23) 

as the first approximation. 

A comment is in order. In place of (123]) it may be also reasonable to take 

£ tH ^ e ^t e tRL +e t{-2iuL3-(ji+v)K 3 +vK + + t MK-} e tKL- 

However, we don't consider this approximation in the paper. 

By the way, we have calculated the term e t{-{^+v)K 3 +uK + + t iK-} - in pj resu ^ j s 

e t{-(ji+u)K 3 +uK++iJ,K-} _ e G(t)K+ e -2log(F(t))K 3e E(t)K. ^ ^4) 

or more explicitly 



RHS =— - exp (G{t)a) ® a T ) (exp (- \og{F{t))N) ® exp (- \og(F(t))N)' 
b (t) V 



x 



exp (E(t)a ® (a f ) T ) (25) 

with 

sinh i^t) sinh (^t) 



cosh (£?t) + ^ sinh (^t) ' cosh (^t) + ^ sinh (^f) 

F(f) = cosh h^-tj + ^ sinh f^tj . (26) 
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On the other hand, we can calculate the term z^l 3 +kL + +k,L-} eas jiy because of the 
relation [L + , L_] = in (TH?1) . The result is 

e t{-2iuL3+RL++KL-} _ e f(t)L +e g{t)L Se l(t)L- = e f(t)(J+-K+) e <?(*)£ 3 e l(t)(J--K-) ^7) 

with 

g— 2iu>t -y Q—liujt j 

/(*) = _ 2 ^ «, g(t) = -2iut, l{t)= _^ k. (28) 

More explicitly 

RHS =exp (/(t) a t (g) (at) 2 ") (exp (-M( a t) 2 ) ® exp f-M(( a t) 2 f 

#0) ,A ^ „„„ ( 9(t) -\jT 



exp ( ^y~N \ ®exp I -^iiV ] x 



exp (l(t)a g> a T ) (exp (~^a 2 ) ® exp (-^(a 2 ) T ) J (29) 



because of ( 12011 . 

Therefore our approximate solution is 



e ± V 1 *e^"^ + ^ K3+l/K + + ' tK -^e ti " 2 ^' L3+Sij++K ' L - | p(0) 



t„t{-(v+v)K 3 +vK + +pK-} t{-2iu)L 3 +iiL + +K.L-}' 



e 2 



exp (G(i)a t ® a T ) x 



F(t) 

(exp (- \og(F(t))N) <g> exp (- log(F(t))A^) T ) x 
exp (F(t)a® (a f ) T ) x 

exp (/(t)at ® ( a tf ) (exp (-^(a^) ® exp f-M ((a t 



\2\T 



2 v y J "V 2 
exp (^^^J ®exp I -^N T J x 

exp (l(t)a® a T ) (exp (-^a 2 ) ® exp (-^(a 2 ) 2 )) p(0) 

and we restore this form to the usual one by use of The result is 

= e -^tjJ2^T (« t ) n {exp({-log(F(t))}iV) x 



n=0 

£ ^ : a m 0W(« t ) m |> exp ({- log(F(t))}iV)K (30) 



x 



,771=0 
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and 




(31) 



This is indeed complicated. 

To construct the general solution to the equation (j3J) is very important in not only Physics 
but also Mathematics. However, it is not easy at the moment, so we only constructed some 
approximate solution. In the very near future we would like to do it. 

In this paper we revisited the quantum damped harmonic oscillator with generalized Lind- 
blad form and constructed some approximate solution in the operator algebra level. 

The model is very important to understand several phenomena related to quantum open 
systems, so the general solution is required. 



On the other hand we are studying some related topics from a different point of view, see 
[8] and [9]. 



Lastly, we conclude the paper by stating our motivation. We are studying a model of quan- 
tum computation (computer) based on Cavity QED (see [TOj and [H]), so in order to construct 
a more realistic model of (robust) quantum computer we have to study severe problems coming 
from decoherence. 

For example, we have to study the quantum damped Jaynes-Cummings model (in our 
terminology) whose phenomenological master equation for the density operator is given by 







i[Hjc, p] (a)ap + pa) a — 2apa)^ (aa) p + pact) — 2a) pa) 



(32) 
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where Hjc is the well-known Jaynes-Cummings Hamiltonian given by 



Hjc = "7T a 3 ® 1 + u o^2 ® c^o + ^ (o+ ® a + cr„ ® a 1 ") (33) 



2 

flat 




with 




\ / 1 \ / 1 

0+ = I I , = I ' a 3 = , 1 2 = 

loy \ o -l y V 01 

Note that p e M(2; C) ® M(.F) = M(2; M(^)), where M(.F) is the set of all operators on the 
Fock space T. See for example [12], [13]. 

Furthermore, it may be possible to treat the generalized master equation given by 

—p = —i[Hj C , p] — 7^ (e^ap + pa} a — 2apa^ — — (aa) p + paa) — 2a 1 pa) 

~ (a 2 p + pa 2 - 2apa) - ~ ((a f ) 2 p + p(a f ) 2 - 2a f pa f ) (34) 

with the condition p,v > \k\ 2 similarly in this paper. 

These equations ((1321). ( 1341) ) are very hard to solve in the operator algebra level, so even 
constructing approximate solutions is not easy. This is our future task, 
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